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Fis  I.  Discs  Posessing  Rotary 
But  No  Approach  Motion  With 
Entrapped  Fluid  Weooe. 


THE  TRANSIENT  LOAD  SUPPORTING 
CAPACITY  OF  FLUID  FILMS 

by 
&?"£.  Gatcombe  (1) 
A  study  Is  presented  In  this 
paper  of  the  thickness  and  the  load 
supporting  capacity  of  fluid  films 
entrapped  between  rotating  circular 
discs,  under  transient  conditions. 
Part  one, (see  Fig.  l)  is  devoted  to 
a  theoretical  analysis  of  the  thick- 
ness and  the  load  supporting  capacity 
of  fluid  films  under  steady  state 

conditions.   Part  two  presents  the  results  of  an  experimental 
investigation  of  these  thicknesses  and  loads  under  transient 
conditions. 

Car-ful  instrumentation  of  this  very  complex  non-steady- 
state  problem,  has  yielded  solutions  which  have  not  been 
obtainable  by  theoretical  studies,  except  for  cases  involving 
assumptions  which  greatly  alter  the  true  results.   It  is  found 
that  the  approach  motion  of  the  surfaces  bounding  the  film  is 
responsible  for  a  much  greater  part  of  its  load  supporting 

(l)  Associate  Professor  of  Naval  and  Marine  Engine extlng,  The 
United  States  Naval  Postgraduate  School,  Annapolis,  Md, 

From  a  portion  of  a  research  project  on  the  study  of  the 
load  supporting  action  of  fluid  films  sponsored  by  the 
Bureau  of  Ordnance,  Washington,  D.  C. 
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action  than  is  the  steady  state  rotational  motion. 
STATEMENT  OF  PROBLEMS 

The  main  problem  is  the  investigation  of  the  thickness 
and  the  corresponding  load  supporting  capacity  of  fluid  films 
entrapped  between  circular  rotating  discs  on  parallel  shafts 
which  have  their  center-to-center  distance  changed  so  as  to 
cause  the  film  to  become  alternately  thin  tl    hick, 
METHOD  OF  ATTACKS 

An  exact  analytical  solution  to  the  steady  state  problem 
was  secured  (see  Part  one).  The  author  had  previously  obtained 
an  approximate  solution  (2), 

The  transient  state  problem  was  solved  by  experimental 
means;   it  is  much  too  comolicated  to  solve  by  analytical 
methods  if  one  retains  such  factors  as  the  deformation  of  tl 
surfaces  of  the  discs,  surface  roughness, .etc.  These  factoi 
occur  in  their  natural  way  in  the  experimental  solution. 
DESCRIPTION  OF  THE  EXPERIMENTAL  SET-UP  (3)s 

Two  parallel  steel  shafts  (see  Fig.  2E)  were  pre; 
and  mounted  with  adjustable  taper  roller  bearings  in  hangars, 
on*  of  which  was  attached  to  the  upper  beam  with  a  strong-bs 
with  an  adjusting  screw  for  controlling  the  initial  gap 
between  discs.  The  other  shaft  was  mounted  in' adjustable 

roller  bearings  in  hangars  secured  to  the  lower  beam  which 

(2)  See  Bibliography 

(3)  See  Fig,  IE 


.could  be  vibrated  transversely  by  means  of  a  Westinghouse  Fatigue 
motor.  The  frequency  and  amplitude  of  vibration  were  controlled 
by  an  oscillator  and  an  amplifier.  No  lost  motion  was 
permitted  anywhere.   Immediate  detection  of  lost  motion  was 
aveilable.  Each  shaft  was  driven  by  8  separate  air  turbine 
wheel  mounted  on  the  shaft.  A  D.C.  excited  coil  was  placed 
in  the  space  provided  on  the  lower  shaft.  A  pick-up  coil  was 
arranged  on  the  upper  shaft.  These  coils  did  not  revolve  with 
the  shafts.  As  the  lower  beam  vibrated  a  signal  occurred  in 
the  pick-up  coil.  This  signal  was  fed  directly  to  a  galvanometer 
type  oscillograph  fitted  with  an  automatic  camera  which  gave 
permanent  records  of  the  signal  generated.  Oscillograms 
were  taken  for  cases  of  the  discs  revolving  with  no  oil  on 
their  surfaces,  a  certain  height  of  wave  form  on  the  oscillo- 
gram results.  Then,  another  oscillogram  was  secured  without 
changing  a  setting  of  any  kind,  but  permitting  oil  to  flow 
onto  the  surface  of  the  lower  discs  on  the  oncoming  side. 
The  second  oscillogram  showed  a  much  reduced  amplitude  of 
vibration  of  the  be8m«  The  thickness  of  the  film  and  its 
load  could  then  be  determined  by  methods  outlined  later. 
Several  sets  of  such  records  were  taken  for  different  settir 
of  the  amplifier. 

Lead  wires  were  rolled  between  the  discs  a  sufficient 
number  of  times  to  insure  no  "spring  back"  of  the  thickness 
of  the  leads  with  the  lower  beam  supported  by  jacks  to 


4 
prevent  bending.  The  thickness  of  these  leads,  determined  by 
interferometer  methods,  gave  the  neutral  gap  thickness  h0„ 
Other  leads  were  rolled  between  the  discs  and  pounded  by  the 
vibrator  until  the  height  of  the  scope  wave  form  reached 
that  which  it  had  when  the  gap  contained  sir.   These  lead 
thicknesses  gave  the  amplitude  of  vibration  of  the  lower 
shaft  discs  for  various  power  inputs  to  the  vibration  n 
coils  and  the  pick-up  excitation-coil. 

The  pick-up  apparatus  showed  a  nearly  linear  response 
for  these  small  amplitudes  of  vibration -(0,001  approx.). 
The  dynamic  spring  constant  of  the  beam  was  secured, 
there  resulted  a  means  of  obtaining  a  continuous  picture  of 
the  thickness  of  the  film  and  its  load  supporting  capacity 
under  these  transient,  conditions. 

The  existing  experimental  set-up  could  not  be  used 
to  obtain  a  steady  state  solution  of  the  problem,  since  the 
pick-up  apparatus  did  not  respond  well  to  A.O.  excitation. 
A  relatively  large  neutral  gap  setting  was  chosen9  which, 
according  to  existing  theory,  would  yield  almost  no  steady 
state  load  supporting  component.  Since  there  is  little 
state  load  for  the  particular  chosen  neutral  fap  settin 
beam  may  be  considered  as  vibrating  about  its  neutral 
position. 


nomenclature: 

a  =  pole  distance 

a     s  abscissa  of  the  bl-polar  coordinate  system,  pressure 
angle 

P •■  ■  ordinate  of  the  bi-polar  coordinate  system 

f  ■  stream  function 

d  ■  distance  from  the  origin  to  the  center  of  the  disc  in 
the  bi-polar  solution 


0  .  12^a  P  V2rho 

h*  «  reciprocal  of  "stretch  factor" 

&Q  ■  neutral  gap  thickness  between  discs 

p  *  pressure 

P  m  load- supporting  capacity  of  the  film  per  unit  face  width, 
pitch  point 

r  s  rlr2 


r  ♦  p 
1    r2 


r^  and  r2  -  radii  of  curvature  of  the  discs 

8  *  tan"1      x   

V  2  rho 

T  *  torque 

U^  and  U2  fc  peripheral  velocities  of  rotating  discs 

U  «  Oj  +  u2 
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u,  v,  and  w  ■  coordinate  velocity  components 


Xq  «  effective  width  of  band  of  contact 

-2    $2  $2  d2     =  Laplscian  operator 

V  *  dx2     dy2     d*2 

jia  *  viscosity  at  atmospheric  pressure  and  stated  temperature 

yx       «  viscosity 

CO      «  angular  velocity 

p      -     density 

a   __  du   +   dv   +    dw 
e   *  "Si   *   3y        Si 

JX  *  hf   by  definition 

INTRODUCTION; 

For  some  time  the  author  has  been  aware  of  the  fact  that 
the  transient  state  of  affsirs,  in  many  hydrodynamical  problems, 
accounts  for  the  relatively  large  load  supporting  capacity  of 
fluid  films.  «ost  of  our  practical  engineering  applications, 
involving  fluid  films  such  as  the  lubrication  of  various 
bearings  end  gear  teeth,  come  under  this  category.  However, 
these  particular  non-ste^dy-state  problems  are, at  present,  wholly 
unsolvble  by  analytical  means,  at  least  in  their  entirety. 

It  has  been  only  recently  that  we  at  the  Naval  Post- 
graduate School  have  acquired  accurate  and  delicate  enough 
equipment  to  undertake  the  solution  of  such  a  complex  problem. 
fart  two  of  this  report  presents  the  results  of  en  investi- 
gation of  the  load  supporting  capacity  of  two  oils,  namely, 
Mavy  §160  and  2190T.  The  results  have  been  limited  by  the 
capacity  of  our  equipment,  but  they  show  clearly  thpt 


relatively  thin  oil  film  loads  and  thicknesses  can  and  Viave 
been  measured  by  these  methods.  Modifications  are  already 
under  way  which  will  permit  the  investigation  of  heavier 
film  lords. 

We  now  focus  our  attention  on  the  analyticsl  solution 
of  the  steady  state  problem.   The  bi-poler  coordinate  system 
is  conveniently  employed,  since  it  provides  boundaries  of 
exactly  the  nature  required  in  this  investigation. 
PROCEDURE: 

The  general  procedure  followed  in  the  solution  of 

this  problem  is  to  {a)    write  the  equations  of  viscous 

fluid  motion  in  the  oartesian  coordinate  system,  (b) , 

assume  that  the  fluid  is  incompressible  and  that  constant- 

viscosity-steady-state-conditlons  exist,  (c)  accomplish  the 

transformation  to  the  bi-polar  coordinate  system,  id)  obtain 

the  expressions  which  indicate  the  characteristics  of  the 

fluid  film. 

EQUATIONS  OF  MOTION  FOR  STEADY  STATE  CONSTANT  VISCOSITY 
INCOMPRESSIBLE  FLUID  CONDITIONS 

Fig.  2  shows  an  elementPl  cube  of  t^e  fluid  with  its 


indicated  positive  stresses.   Texts 
on  hydrodynamics  present  the  follow- 
ing equations  for  viscous  fluid 
motion  in  which  the  viscosity  p.   is 
considered  constant. 
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If  the  body  f trees  X,  Y,  and  Z,  are  neglected,  the  fluid 
assumed  Incompressible,  and  steady  state  exists,  then 
Equations  (l)  reduce  to 
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If  the  equations  of   (2)  are  differentiated  with  respect  to 
x,  y,  and  z  resoectively,  and  added,  then  the  following 
equation  results! 
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We  ere  now  concerned  with  two  dimensional,  flow  c  iy,«o 
Equation   (3)  may  be  expressed  aa 
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Equation  (3)  la  now  transformed 
into  an  equation  in  bi-polar 
coordinates.   As  has  been 
previously  mentioned,  these 
coordinates  permit  boundary 
conditions  which  are  exceedingly 
suitable  to  this  problem.  Ho 
approximations,  in  regard  to 
film  thicknesses,  need  be  made. 

The  equation  which  is  used  to 
transform  from  C8rtesian  coordinates  _ 
to  bi-polar  coordinates  is 


r;-ac 


Fig. 3  Circular  Discs 
Associated  With  The 
Bi-Polar  Coordinate 
System 


(5) 


Pig.  3  presents  the  ac,  p  plane  of  the  bi-polar  coordinate 
system  superimposed  upon  the  x,  y,  plane  of  the  Cartesian 
system.  The  origin  for  both  systems  is  at  0.  The  equation 
of  a  circle  which  represents  the  circular  disc  of  this 
problem  is  as  follows: 


X  +  (y-   a  co/ac*^    af'csch  <K 


(6) 


The  distance  (d)  from  the  origin  to  the  center  of  the  circle 
is  given  by 

d  =  a  coth     <L     (7) 
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Smell  a  represents  the  distance  from  the  origin  to  the  focal 
point  of  the  bl~polar  system.   The  radius  (r^)  of  the  disc  is 
expressed  by 

r*i      a  c sen  cl   •  ••••••••e««*.*»«*«..«. ...•*•  \ o / 

The  fc   owing  relationships  exist* 

a  sin  n.h  a 

X   =  — — — —  jy   s  _— . „..(9) 

cosh  «c  "  -os  (3  cosh      os  p 

.ditioi   of  an  Incompressible  fluid  we  may  write 

du        dv 

o-        dy 

We  may  also  write 

^V         3^        S>*-   a„ 

M  dx  ^   J        ^      s^       rT       l  y 

where  I  is  the  stream  function. 

In  the  Cartesian  system  we  may  write 

r-Fr-    -aJ:^  -fer— "       - ^ 

In  the  bi-polar  system  Equation  (12) Is  written  as 
A. 


where  (h)   is  the  reciprocal  of  the  'stretch  factor™. 
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It  Is  expressed  symbolically  by 


_  -     L-_\^-  / )  _. _ ,(if) 


or 

cosh  a  -  cos  p 


(15) 


2 

The  differential  equation  satisfied  by  ^ *  is 

^         ^     -  O  (16) 


,,,^ 


An  appropriate  solution  is  (4) 


where  the  coefficients  An  ,  Bn  ,  Cn  ,  %  are  arbitrary 

constants,   if  we  write  h¥  =Y  ,  the  corresponding  expression 

for  ?£  may  be  determined  with  the  aid  of  the  differential 

equation 

(4)  See  Bibliography 


CL%=    (C.sk«-  C+s/3)(  ^  ^  ^ 


v- 


-+-  (cosh  <<-f-&>s/5)X* 

It  is  found  thUt 

2  CasVi/3//)v,0^(H^«H  +&>>  Cosh  (*-/)*{ 

Sine©  the  quantities?  end  -//J*  are  conjugate  function, 
equations  of  the  following  type  are  to  be  satisfied} 

2z    jif  ^y    ~ - — o&) 

Thus  an  expression  for  p  in  the  a,  p*  system  is 

Equation  121)  represents  the  general  expression  for  the 
pressure  p.   If  we  consider  problems  in  which  the  press 
is   single  valued^then  all  terms  possessing  d   or  p  alon< 
will  not  be  retained.  The  coefficient  (B0  +  I>i)^of  p  may 
set  equal  to  zero. 


13 


Brt   +  D,    a   0 


(22) 


If  Utf  ,  up   ere  the  velocity  components  in  the  positive  cr, 
and  p  directions,  we  have 


and 


u«  -  -  *  If  -  -  f  + 1    •*  p 


UB  =   ~  h  35  =   oa  "  a   flinh 


(23) 


(24) 


We  let  the  angular  velocities  of  the  cylinders  be  represented 
by  t*» and  write  the  following  equations  which  represent  the 
boundary  conditions* 

dX-    Yi 

For  d  =  +di,   T  =  Yi,    and  a£  =  a      sinh  ^i   "   a^eosech  *1      ....(25) 

For  ci   =  -ax, *  =  -*i,  and  |£*  =  -J3-  sinh  (-cti)  -  a^cosechl-ai) 

(26) 


<3X* 


By  symmetry  for  ct  =  0 ,   Y  =  0„   and  -r—g   = 


0 


(27) 


Thus  the  expression  for  }Cof  Equation  (19)  may  be  reduced,  in 
the  case  of  our  problem,  to 


/C  =  B^cC(cosh  <£  -  cosS)  +  C  sinbjc  +  C  sinh2et  cos{3 


(28) 


where 


Z  Cos  A  <*,   5tnh\t 

dw 


Cr- 


(29) 
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If  this  simplified  expression  f or)£ ,   Equation  K 28)^  is  substituted 
into  Equation  (18).  the  following  expression,  appropriate  to 
this  problem,  may  be  written* 


*~4C,(* 


— -C3e) 


C.  foink*  c*>s0-2,$mh^C*sl\*<  CoS/5  f-$tnkc(  Cosh<* 

The  derive tive  with  respect  to  a   yields 

&  §^-  -f  C,  [cos/iflf  CoT^  -  *  Cotp  Cosh  2o(i-  Cosh  Z<<  )-—  (?0 

Integrating  partially  with  respect  to  (3  following  the  statement 
of  Equation  ( 20) gives 

f 

CL-ff  «  C  +  '^/*  C,  Co*ii*  Sm£ •  £ c*shz«  S*» 2/S    —--  pz) 

Mow  a  comb 5   tion  of  facts, t  hst  there  is  symmetry,  that  the 
fluid  is  assumed  in^spable  of  exerting  tension  stresses,  that 
we  are  dealing  with  conversing  and  diverging  spaces  between 
the  discs,  require  thst  the  pressure  shall  pass  through  zero 
at  the  oriscin  where  (3  =  ir. 

This  yields  G  -  0.      Thus  our  final  expression  for  the  pressure 
p,  in  the^l  (3  system,  is 

-y>*  &£££costioC3t*0  -  £  Cosk  JL*i  Ss*<  X&J- 

where 


pO 


a.  *+> 


4c»*4'<l3'#ti*>(i 
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An  expression  for  the  maximum  pre a aura  (Png^  )»  along  the 
x-axis  or  a  *  0,  It 

^    m      3)/T^/U  C, (36) 

an  expression  is  now  obtained  for  the  load-carrying 
capacity  (F)  per  unit  face  width  of  the  fluid  wedge.  We  note  from 
Fig,  3.  that  an  expression  for  this  force  is 

<*°  -  /*  f  c*s  4  c/6  (58J 

O 

It  is  necessary  to  express  p  of  Equation  36  in  terms  of  #, 
remembering  that  «r  may  be  considered  constant  for  any  given 
sice  disc.  Sxpressions  for  sinp  and  ain2p  are  the  following: 


3t«  s#    -    ?*'"  4*"*  **/     (37) 

cL-t-r  Cos  <p 

( cj  -t  r  cos  ff~ 
in  expression  for  F  is 


? •   Z^C' %"*&'  (*/  * *z)  (39) 


where 


Kf=  ^(2^^^-^W^^, 


(40) 


and  K=2qs,»£x<:*'4x-<,Uts^-  2tc/   7      (41) 

We  further  observe  that  nearly  all  the  load,  carried  by  the 
fluid  wedge,  ie  distributed  over  a  very  small  area.  Let  us 
arbitrarily  define  the  width  of  the  supporting  band  of  contact 
xQ   as  that  distance  from  the  origin,  measured  along  the  x-axis, 
to  the  point  where  the  pressure  drops  off  to  one  per  cent  of 
its  maximum  value  *  An  approximate  expression  for  x0  is 

x0   ■  6.3  a . . «.(42) 

From  previous  work  by  this  author  (5)  it  was  found  tha 
an  approximate  solution  to  this  problem  gave  the  following 
set  of  equations* 
For  the  pressure  £ 


1"  g(- 


.h.,.  /:„  '**"Y77hT .(44) 


U  is  the  peripheral  velocity  of  the  disc. 

~JL^ — — -       For  equal   sise  discs  r  »  -*»»       ...... 

*1      *       *2 


*s~V2rhl    tan.  S 

o 

(5)  See  Bibliography 


•  .  .  s  a>  44©  / 


17 

For  the  maximum  pressure  P^a* 

-ft".*  -7F  G—.r-r - -  (+V 

For  the  load-carrying  capacity  P 

f>=.  ^VZTJu    (+8) 

For  the  width  of  the  band  of  contect 

x.-YzrXrio'iC8S*sy~ (4ff 

•e  shall  now  apply  these  equations  In  analyzing  a 
given  problem.   Since  the  exact  solution  is  obtained  for  the 
case  of  symmetrical  flow. only,  we  shall  take  discs  of  equal 
size  revolving  at  the  same  speed  CV      „  The  discs  are 
3.0CC0  inches  in  diameter.   Their  sneular  velocity  is  (-183) 
red/ssc.   The  minimum  thickness  of  the  wedge  Is  h0   =  0.00002*BI 
Inches;   The  viscosity  of  the  oil  lsjL/    ■  5.5  x  10"6  lb   sec/sq.ln 
We  shall  compute  the  following  by  the  eatsct  and  approximate 
methods  and  compare  the  results: 

1.  The  maximum  pressure  within  the  wedge 

2.  The  load-c8rrylng  capacity  P  of  the  wedge  per  ino.h 
of  face  width 

3.  The  width  x0  of  the  supporting  bsnd  of  contact 

SOLUTION  OF  PROBLEM;  -  EXACT  METHOD 
Equation  (35)  is  the  expression  for  the  maximum' pressure 

iW- sr— lSo> 

_fi_  =  -  -i S* -, — -(St) 


18 


To  find  ai   we  need  to  use  the  following  z 
d-^  =  a  coth  GLi   ;  f?     =  a  csch  a^ 

or  &     =  cosh  c*i 

d  =  1.5000   +  0.00001«*  =   1.50001«i 
so   cosh  ax     »  1.000006667 
c£l     =  0.00365  red. 
and  a  s/'sin/ic^     =  1.5000  x  0.00365 
a  *  0.00547 


inches 


Thus 


15^*-  sYT  xS'Sx/* 


-e 


'MA** 


/03 


#X/'  »e«®  GG€7X  O.  &o  9CS 


-...(SZ) 


%4*.  "  +  *?C*°  ■#•«*' 


Equation  (39)  is  the  expression  for  the  load-carrying 
capacity  F 


Kx  =  +0.00084 
Kg  =   -5.97 

F  =  -226  lb  in.  of  face  width (55) 

Equation  (42)  is  an  expression  for  the  width  of  the  supporting 
band  of  contact 
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x0  =  6.3a  =  6.3  x  0.00547  ■  0.0344   Inches .(56) 

SOLUTION  OF  PROBLEMS    -   APPROXIMATE  METHOD 


The  maximum  pressure  P^ax     is   given  by 


1*      *   -t-^-  G    —  -Z-Z--XZ+7  +  .*  ~  2GM  1*o  7*>* 


t 


« 


It 


The  load-carrying  capacity  £  is  indicated  by 

JP      =  225  lbs.  per  inch  of  fece  width. 
The  width  of  the  bend  of  contact  is 

COMCLUSIOH  COHCSRHINQ  STEADY  STATE  SOLUTION 

We  find  that  the  approximate  solution  derived  earlier,  gives 
results  which  are  close  to  the  exact  answers.  However,  at  this 
point ,  our  results  are  of  little  account  in  the  solution  of  the 
usual  non<=steady-state  problem  except  for  a  guide  in  the  solution 
of  steady  state  problems.  If  we  wish  to  get  practical  results 
for  these  problems  In  hydrodynamics  we  are  forced  to  get  the 
solution  of  these  non- steady- state  problems  through  experimental 
analytic  We  now  proceed  to  Part  two  of  this  report. 


PART   II 
INTRODUCTION; 

The  thickness  and  the  load  supporting  cape city  of  oil 
films  entrapped  between  circular  rotating  discs  on  parallel  shafts 
were  investigated  here  by  experimental  methods.   (see  Fig.  IE). 
The  fluid  wedge  was  formed  by  dragging  the  oil  between  the 
discs  which  were  rotated  in  opposite  directions.   The  wedge  was 
made  alternately  thick  and  thin  by  vibrating  the  beam  which 
supported  the  lower  shaft.   The  D.C.  excitation  coil  set  up 
lines  of  flux  which  had  to  pass  through  the  neutral  gap  which 
separated  the  surfaces  of  the  two  sets  of  double  discs.   The 
pick-up  coil  delivered  a  signal  to  the  galvanometer  oscillograph 
as  the  lower  shaft  was  vibrated.   An  automatic  camera  was  used 
to  photograph  the  oscillograms.   Oscillograms  were  taken  with 
and  without  fluid  wddges  between  the  discs.   The  presence  of 
oil  caused  the  beam  supporting  the  lower  shaft  to  vibrate 
through  rel"  ively  small  amplitudes.   Small  lead  wires  of  the 
order  of  0.020  inches  were  rolled  between  the  two  sets  of  discs 
with  the  lower  beam  supported  by  jacks  to  prevent  bending. 
The  thickness  of  the  neutral  gap  h0  was  measured  by  obtaining 
the  thickne.ss  of  the  flattened  lead  wire.   The  neutral  gap 
thickness  ho  was  adjusted  until  the  amplitude  of  the  vibrating 
beam,  for  a  power  setting  of  the  amplifier  of  60,  equaled 
the  value  of  the  gap  h0. 

The  constant  D.C.  excitation  current  for  the  pick-up 
system  came  from  two  storage  batteries.   The  D  C.  and  A.C. 
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currents  and  voltages  to  the  fatigue  motor  were  fixed  for 
various  amplifier  settings.   Then  the  values  of  the  height 
of  the  wave  form  on  the  scope  were  found  for  various  settings 
of  the  amplifier  with  the  angular  velocities  of  the  discs  at 
zero.   These  values  are  plotted  in  tfraph  No.  1.   Lead  wires 
were  then  rolled  between  the  discs,  left  in  position  with 
co  =  0,  and  the  vibrator  pounded  them  until  the  height  of 
the  scope  wave  form  equaled  that  which  existed  for  the  case 
of  air  between  the  discs.  Leads  were  taken  and  their 
thickness  measured  for  all  six  power  settings  of  the  amplifier, 
except  no.  50  where  the  leads  became  too  thin  to  handle. 


EXPERIMENTAL  RESULTS: 

It  is  clear  now  that  at  this  point  we  have  Information 
at  hand  for  determining  the  magnification  factor  (M.F.)  of 
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our  set-up.   For  example,  we  find  from  Graph  No.  1,  that  for 
an  amplifier  setting  of  60  (without  oil  and  co  =  0)  the 
lower  disc  vibrated  through  an  amplitude  of  h0  =  0.0C197 
inches,  while  the  height  of  the  scope  wave  form  v  as  1.80 
inches.   This  s-ives  a  M.F.  =  914.   Likewise,  when  the  amplifier 
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was  set  at  30  the  leads  showed  that  the  lower  shaft  vibrated 
through  aa  =  0.0C099  inches,  while  the  scope  gave  as  = 
1.44  Inches.   The  M.F.  for  this  setting  was  1450. 

It  was  noticed  that  the  as  depended  not  only  upon  the 
amplitude  of  vibration  of  the  disc  a^  ,  but  also  upon  the 

speed  with  which  the  discs  rotated.   It  thus  became  necessary 
to  plot  the  scope  amplitude  against  the  amplifier  settings 
for  various  values  of  the  angular  velocity  Gc?  .   These 
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25 
value a  are  shown  In  Graph  No.  1.  The  magnification  factor 
M.F.  will  always  be  determined  with  reference  to  the  graph 
of  co  =   0.  The  increase  in  the  amplitude  of  the  scope  wave 
due  to  the  rotation  of  the  discs  is  nearly  nullified  for  the 
cases  with  and  without  oil,  since  the  angular  velocity  was 
fixed. 

We  indicate  now  how  one  may  use  Graph  No.  1  to  determine 
the  load  carried  by  a  certain  thickness  oil  film  wedge.  Prom 
Table  No.  1  we  find  the  following:  that  the  height  of  the 
scope  wave  form,  corresponding  to  w  =  2400  RPM  without  oil 
with  the  amplifier  at  60,  is  2.08  inches,  that  the  amplitude 
of  the  scope  wave  form,  corresponding  to  «*»  ■  2400  RPM  with 
N3100  oil  with  the  amplifier  at  60,  is  1.35  Inches.   We  note 
that  0.28  inches  of  the  scope  amplitude  are  due  to  rotational 
effects.    ->w  the  amplitude  of  vibration  of  the  lower  disc, 
without  oil,  with  the  amplifier  at  60  and  f  =  425  cps  is 
just  sufficient  to  close  the  neutral  gap  setting  of  hQ  ■ 
0.00197  inches.   Thus,  for  the  case  of  no  oil,  with  w   =   2400 
RW,  with  f  -  425  cps,  the  M.F.  is  found  to  be  1.80/ 0.00197 
or  914.  The  magnification  factor  is  defined  as  the  ratio  of  a  / 
ad  .   This  value  is  calculated  with  reference  to  the  curve  ^ 
=  0.  For  the  case  of  the  use  of  N3100  oil  with  <*    «  2400  RPM, 
with  f  =  425  cps  we  find  from  Table  one  that  the  a8  =1.35 
inches.   We  recall  that  0.28  inches  of  this  amplitude  is  due 
to  rotational  effects,  thuB  the  value  corrected  to  the  <*»  «  O 
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TABLE     No. 

1 

Nov.  1950 

RPM 

Temp.        70°  F. 

No  Oil 

2I90T 

N  3100 

NftQjl 

AMPLIFIER 
SETTINGS 

as  ins. 

osins. 

ad  ins. 

M.F. 

!    INS. 

F0  lb. 

Fs.s.  lb. 

as  ins. 

ad  ins. 

M.F. 

t     INS. 

F0  lb. 

F..  lb. 

P0  lb. 

60 

2.0  8 

1.94 

0.00127 

1360 

0.00070 

140 

20 

1.35 

0.00081 

914 

0.00116 

162 

73.3 

102.6 

50 

2.04 

1.82 

0.001 1  1 

1420 

0.00086 

140 

1.25 

0.00073 

1280 

0.00124 

158 

98.4 

40 

1.88 

1.80 

O.OOMO 

1420 

000087 

140 

1.17 

000063 

1268 

0.00134 

156 

93.6 

30 

1.70 

1.50 

0.00087 

1350 

0.00110 

1  17 

0.97 

0.00061 

1  156 

000136 

134 

81.9 

20 

1.05 

088 

0.00054 

1  146 

000143 

69 

0.63 

0.00045 

976 

0  00152 

74 

48.9 

10 

0.45 

0.39 

0.00043 

635 

0.00154 

21 

0.36 

0.00035 

622 

0.00162 

22 

20.4 

<*>    =  1750    RPM 

No  Oil 

2190  T 

N3I00 

No  Oil 

AMPLIFIER 
SETTINGS 

a,  ins. 

a,  ins. 

ad  ins. 

M.F 

t    INS. 

F0  lb. 

F,.,.lb. 

a,  ins. 

ad  ins. 

M.F. 

t      INS. 

F0  lb. 

Fl,.,,b 

Pf  lb. 

60 

2.00 

1.82 

000126 

1341 

0  00071 

136 

15 

114 

000072 

1255 

0.00125 

164 

53.3 

102.6 

50 

1.90 

1.69 

000100 

1445 

000097 

142 

1.12 

0.00069 

1252 

0.00128 

164 

98  4 

40 

1  85 

1.65 

0.00099 

1450 

000098 

132 

105 

0.00064 

1250 

0  00133 

154 

93  6 

30 

173 

1.30 

0.00076 

1309 

0.00120 

123 

0.89 

0.0005  8 

1214 

0.00139 

133 

81.9 

20 

1.00 

0.75 

000057 

1060 

000140 

68 

0  59 

0.00046 

976 

000151 

73 

48.9 

10 

0.40 

0  30 

0.00043 

620 

0.00154 

20 

0.29 

000039 

564 

0.00158 

20 

20.4 

^    =  875    RPM 

No  Oil 

2190  T 

N  3100 

No  Oil 

AMPLIFIER 
SETTINGS 

a,  ins 

osins 

ad  ins. 

M.F. 

t    INS. 

F0  lb. 

FSJlb. 

a,  ins 

ad  ins. 

M.F 

t    INS 

F9lb 

*Y»'b 

P9lb 

60 

1.80 

1.74 

0.00129 

1336 

000068 

136 

7.3 

1.07 

000079 

1250 

0.001  18 

164 

26  7 

102.6 

50 

1.78 

1  68 

0.001  1  8 

1405 

000079 

134 

1.04 

000077 

1320 

0.00120 

154 

984 

40 

172 

1.65 

0  001  17 

1404 

000080 

133 

0  95 

000070 

1262 

000127 

150 

93.6 

30 

1.44 

1.40 

0  001  10 

1320 

0.00087 

105 

0  82 

0  00064 

12  50 

000133 

154 

819 

20 

080 

0.84 

000064 

1250 

000133 

65 

0  55 

000047 

1046 

000150 

73 

48.9 

10 

0  33 

0.30 

0  00043 

620 

0.00154 

20 

0.26 

0.00037 

560 

000160 

20 

20.4 

DYNAMIC    SP    CONSTANT  OF   BEAM  =    250,000    lbs  /IN. 

VOLTAGE    TO  VIBRATION    MOTOR     =    8V   (ACCORDING   TO   LOAD). 

CURRENT  TO  VIBRATION    MOTOR    =  200   MA.  (ACCORDING  TO  LOAD). 

VOLTAGE   TO  VIBRATION    MOTOR     =    118  V  DC. 

CURRENT  TO   VIBRATION    MOTOR  =  0  9  A.  DC 

VOLTAGE    TO   PICK-UP    SYSTEM     =    12V.    DC. 

CURRENT  TO   PICK-UP  SYSTEM    =  0  25  A   DC 

t=he-Od    MINIMUM    OIL   WEDGE    THICKNESS. 

TOTAL  WT.  ON   BEAM    =   12.11    lbs 

NATURAL   FREQ    OF   BEAM    WITH  WT     =  12.11    lbs    ■  450  cps 

F0  *  LOAD   CARRYING    CAPACITY   OF   FILM    IN   Ibs/IN.  OF  FACE    WIOTH 

F,,»LOAD   CARRYING    CAPACITY  OF   FILM(FROM    STEADY    STATE    SOLUTION) 

P0  =  AMPLITUDE    OF   MOTOR   OUT   PUT    FORCE. 

No  Oil  f  =  425  cps 
"    •  0    RPM 

.     Ampl   at  60                          s7'fl                                         .'■'■'W 
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27 
curve,  Is  1.07  inches.   This  gives  a  M.F.  =  1326.   This 
scope  amplitude  of  1.07  inches  on  the  c^  =  0  curve  corresponds 
to  an  amplifier  setting  of  24,  and  thus  a  disc  amplitude  of 
0.00081  inches.   The  amplitude  of  vibration  of  the  beam  was 
thus  reduced  from  ho  =  0.00197  inches  without  oil,  to 
0.00081  inches  with  oil.   We  know  now  that  the  minimum  oil 
film  wedge  thickness  was  0.00116  inches.   The  load  on  the 
film  will  be  computed  now. 

From  the  mechanics  of  the  problem  we  observe  that  the 
beam  may  be  regarded  as  massless  with  an  equivalent  weight 
placed  at  the  center  of  the  span.   There  is  also  the  weight 
of  the  lower  shaft,  bearings,  hangars,  motor  drive  coil,  etc. 
*<hich  may  be  considered  concentrated  at  the  midsection  of 
the  beam.   This  total  weight  =  12.11  lbs.   We  see  that  four 
forces  act  on  the  beam  in  the  displaced  position,  namely; 
the  beam  restoring  force  FD  ,  the  force  of  internal  damping 
Fd  ,  the  force        oil  film  wedge  FQ,  and  the  force  of 


the  vibration  motor  Fm. 
The  amplitude  and  the 
phase  relationship  of  these 
forces  are  shown  in  Fig.  3B. 
The  sum  of  the  vertical  com- 
ponents may  be  written  as* 


j4i  X 


<f 


Fie,  3£ 


m*»  y0  +     P0  cos  *     -  F0  -     ky0     *     0 


(57) 
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The  value  of  F0,  the  amplitude  of  the  motor  output  force, 

was  found  for  different  amplifier  settings,  by  placing  SR4 

strain  gages  on  a  specially  prepared  drive  rod.   These  values 

are  Riven  in  Table  No.  2.   It  was  found  that  the  amplitude 

of  the  motor  output  force  P<>  depended  upon  the  amplifier 

setting  for  a  given  frequency.  Thus  the  motor  output  force 

was  evaluated  without  the  use  of  oil  between  the  discs. 

The  phase  angle  #  is  shown  to  be 

tan  $     -         b„„  '  ofc>  ,      , 
ST ii ( 58) 

1   -    (  cu2 /     2 

where  b  is  the  damping  constant  of  the  system  and  is  evaluated 

from  the  die-away  wave  form  shown  in  the  preceding  photograph. 

It  was  found  that  b  =  0.95  lb.  sec.^  in.   The  value  of  the 

critical  damping  constant  bcr    =  ^  4mk    was  found  to 

be  bcp  =  177  lb.  sec./  in.  The  weight  on  the  beam  was 

12.11  lbs.  and  the  dynamic  spring  constant  was  250,000  lbs. /in. 

The  undamped  natural  frequency  of  the  system  was  4jj  n   = 

2830  rad.  /  sec.   It  was  found  that  $  ■  5°  -  10  . 

The  amplitude  of  load  F0  carried  by  the  oil  film  wedge 

may  now  be  found  from  Eq.  57.  The  disc  or  beam  amplitude  of 

Vibrations  a<j  or  y0  was  found  to  be  0.00081".  The  mass  m  on 

the  beam  was  0.0314   lbs.  sec.    ^  The  frequenCY  of  vibration 

in. 

of  the  beam  was  425  cps  or  2670  rad. j  sec.   The  value  of  P0 
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was  102.61  lbs.   Thus  the  oil  wedffe,  the  minimum  thickness  of 

which  was  0.00116  ",  carried  a  load  of 

F0  =  -  250,000  x  0.00081  +  102.6  x  0.006  +  09314  x 

2 

2"570~  x  0.00081      (59) 

or 

F0  -  81,00  lbs. /half  inch  face  width  (60) 

FQ  =  162.0  lbs.  /inch  of  face  width   (61) 

If  one  used  the  results  of  Part  I  (Eo.  56)  to  Dredict  the 

width  of  the  supporting  band  of  contact  x0  ,  then  the  average 

pressure  within  the  film  will  be 

p  en  q 

Pav=  T"      "        0.6  x  0.0344      =  4710  osi.(R2) 
One  may  count  on  the  maximum  pressure  being  approximately 
seven  times  the  average.   This  would  give  a  maximum  pressure 
of 

Pmax    =  33,000  psi   (63) 

It  will  now  be  of  considerable  interest  to  compute  the  load 
which  this  same  m'nimum  thickness  oil  film  wedge  will  support 
on  the  basis  of  steady  state  conditions,   kq .  39  or  Eg.  48 
may  be  used  to  obtain  this  result.   In  Iq .  48 
G  =  12  //  a   U  y/   2rh0 


(64) 


>/a  =   1.29   x  10"4reyns  (^3100   r    70°F) . 

U  =   378"  /  sec.    ( 240C    PPM   -    3"    diam.disc) 
ho=  0.00197 
0  =   8130   osi 
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Thus 

0 


P  =  g    V  2rho 


s  8130 

-g —  x  0.054   ■  73.3  lbs./  in.  of  f ace  w  idth  or 

S     =36.6  lbs./  half  inch  of  face  width. 
This  is  the  load  which  the  oil  wedge  could  be  expected  to 
carry  on  the  basis  of  steady  state  conditions.   It  is 
observed  that  the  non- steady  state  solution  indicates  that 
this  same  wedge  can  carry  81  lbs.  The  non- steady  state 
solution  then  shows  that  the  wedge  can  carry  2.22  times  as 
much  load  as  the  steady  state  solution  indicates.  From  an 
inspection  of  Fig.  4E  it  will  be  observed  further  that  the 
ratio  of  the  largest  value  of  F0  to  the  smallest  value  of 
FS8  is     fw  =  6«  This  means  that  the  load  supporting 
action  of  the  film,  resulting  from  the  approach  motion  of 
the  metal  bcmding  surfaces,  can  be  six  times  as  large  as 
that  resulting  from  the  rotational  effects  of  the  discs. 
In  other  words  the  no;i-steady  state  solution  gives  values 
which  are  six  times  as  large  as  those  of  the  steady  state. 
This  is  a  factor  which  we  may  no  longer  omit  in  our  pre- 
dictions of  the  load  carrying  capacity  of  films.   It  is 
no  wonder  that  we  hsve  not  been  able  to  account  for  the 
large  loads  which  bearings  and  gears  actually  carry.   Our 
theory  is  based  on  steady  state  conditions. 

The  preceding  point  is  the  prominent  part  of  this  whole 
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research  project.   The  oil  film  wedse  in  bearings  and  gears 
supports  loads  because  they  operate  under  non-steady  state 
conditions. 
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Information  on  use  of  strain  gagess 

The  amplitude  of  the  motor  output  force  was  determined 
through  the  use  of  strain  gages.   Two  A-7  strain  gages  were 
used  on  the  flat  section  of  the  drive  rod.   A  third  tempera- 
ture compensating  gage  was  employed.   The  strain  indication 
was  received  on  a  cathode  ray  oscillograph,  through  the 
employment  of  a  Baldwin  Southwark  strain  calibration  and 
amplifier  instrument.   Two  gages  vrere   used  so  that  bending 
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of  the  drive  rod  could  be  detected.   It  was  determined  that 
a  strain  of  120  micro  inches  per  inch  corresponded  to  an 
amplitude  ordinate  height  of  C.7  inches  on  the  cathode  ray 
oscillograph  with  the  gain  as  used.  The  cross-sectional 
area  of  the  steel  drive  rod  was  A  =  0.0100  sq.  ins. 

TABLE  II 


Amplifier 

Amplitude 

Ampl itude 

Id's. 

settings 

Gage  A  in 

inches 

Gage  B  in  inches 

60 

1.99 

2.00 

102.6 

50 

1.90 

1.91 

98.4 

40 

1.80 

1.81 

93.6 

30 

1.58 

1.59 

81.9 

20 

0.95 

0.95 

48.9 

10 

0.35 

0.40 

20.4 

Information  on  oils  used: 

N31QQ  -  navy  symbol? 

specific  gravity  at  60  /  60°P  =  0.903 

A.P.I,  degrees    ■   25.2 

Viscosity  seconds,  S  U  V,  at  100°F  =  1318 

n       n        «     n   130°F  ■  492 

"  »  »  "  210°F  =  90.9 
Viscosity  index  (Dean  and  Davis)  =  72.9 
2190T-  turbine  oil 

specific  gravity  at  60/  60 °F   =  0.8844 
A.P.I,  degrees     ■  28.5 
Viscosity  seconds,  S  U  V,  at  100°F  m     420 

it       it         n     it   130°F   ■   193 

it       n        w     it   2io°F  =   5e 

(Engineering  Experiment  Station  Test  1  BDT031) 
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Additional  Remarks 

During  the  operation  of  the  test  equipment,  it  was 
found  that  an  optional  method  of  measuring  the  value  of  h0 
gave  a  good  check  on  those  values  already  found  through  the 
use  of  leads.   All  these  test  results  are  reproducible. 

In  the  test  prosra  m  to  date  the  value  selected  for 
the  nettral  gap  setting  (h0  =  0.00197')  was  relatively 
large  so  that  the  beam  would,  in  effect,  vibrate  about  its 
static  neutral  equilibrium  position,  rather  than  about  some 
displaced  position.   In  future  work,  if  small  values  of  h0 
are  used,  it  will  be  necessary  to  correct  the  film  thickness 
to  account  for  this  fact.   The  material  in  the  steel  discs 
was  SAE  1045  steel  hardened  to  45  RC.   The  surface  roughness 
was  20  millionths  of  an  inch. 

The  irregular  freauency  components  of  vibration 
present  in  the  various  wave  forms  shown  in  the  photograph 
result  from  various  causes  such  as  non-concentricity  of 
the  discs,  non-uniformity  of  the  rollers  in  the  bearings, 
and  non-uniform  speed  of  the  discs.   These  components  do 
not  alter  appreciably  the  results  of  the  test. 
Conclusions! 

A  prominent  cheractpristic  of  the  load  supporting 

capacity  of  fluid  films  has  been  uncovered  in  this  research  work. 
It  is  definitely  shown  that  thin  fluid  films  strongly 
resist  being  squeezed  out  from  between  metal  surfaces. 
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especially  if  the  bounding  surfaces  are  a  poroacbing  each 
other  at  high  velocities.   In  other  words,  the  non- steady- 
state  character  of  the  motion  of  the  hounding  surfaces  of 
fluid  films  is  responsible  for  their  ability  to  support 
relatively  great  loads. 

Our  steady  state  theory  never  has  been  able  to 
account  satisfactorily  for  the  relatively  large  loads 
which  our  bearings  and  gears  carry  every  day.   This  theory 
indicates  unreasonably  thin  film  thicknesses  for  these 
members  when  used  with  heavy  loads  at  low  speeds;   many 
times  these  thicknesses  are  much  smaller  than  the  surface 
rouffhness  of  the  metal  bounding  surfaces.   Neither  can 

the  theory  account  for  the  loads  which  are  carried 
successfully  in  very  high  soeed  gearing  such  as  in  certain 
turbine  gear  units.   In  such  cases  the  period  of  mating 
is  very  short,  and  the  whole  lubrication  ohenomenon  is 
non-steady-state,  without  question.   It  is  highly  probable 
that  the  prominent  point  uncovered  in  this  work  partly 
explains  why  such  gears,  actually  work. 

Two  external  tyoe  cylindrical  discs,  which  formed 
the  metal  hounding  surfaces  of  the  fluid  film  under  investi- 
gation, were  chosen  in  this  work  in  preference  to  actual 
bearings,  because  of  their  basic  character  in  relation  to 
gears  as  w<=ll  ss  to  bearings.   The  results  of  this  work 
may  be  used  in  studies  of  gears  as  well  as  in  Internal 
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and  external  bearings.   For  one  of  the  cases  examined  In  this 
work,  the  non-steady-state  fisrures  collected  from  ectual 
experimental  tests  Indicated  that  the  film  under  Investi- 
gation carried  a  load  six  times  as  lergfi  as  the  steady  sfste 
solution  would  have  indicated. 

No  attempt  will  be  made  at  present-  to  formulate  a 
non-steady-state  theory  to  explain  this  non-steady-state 
phenomenon.   But  it  will  not  be  difficult  to  do  this  once 
we  have  more  test  data  to  draw  from  in  f ermine;  semi-empirical 
theories.   We  should  consider  this  work  b    simple  beginning 
in  the  investigation  of  a  most  complex  problem. 

It  has  been  the  opinion  of  this  author  for  some 
time  that  theory  alone  is  an  insufficient  tool  for  com- 
pletely explaining  some  of  these  complex  lubrication  phenom- 
ena.  Thus  we  at  the  Naval  Postgraduate  School  have  soent 
the  last  three  years  In  acquiring  and  testing  accurate 
instruments  for  carding  out  this  mission.   It  cannot  be 
overemphasized  that  we  have  simulated  actual  operating 
conditions.   We  have  worked  with  actual  steel  discs  with 
real  lubricants.   These  discs  are  deformable;  they  have 
surface  roughness;   they  have  certain  adhesive  and  cohesive 
properties.   The  lubricants  act  in  their  natural  fashion. 
Above  all  we  have  worked  with  non- steady- state  conditions. 

This  problem  should  be  pursued  vigorously  along 
these  combined  experimental  and  theoretical  lines  until 
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enough  date  is  gathered  so  that  we  may  formulate  the  proper 
semi -empirical  theories  to  give  important  er.uicU.nce  to 
designers  in  the  fields  of  bearing  and  gears  as  well  as  in 
other  related  fields. 
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